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The influence of a uniform rotation with frequency Q. on the critical behavior of liquid *He near T\ 
is investigated. We apply our recently developed approach [R. Haussmann, submitted to Phys. Rev. 
B] which is a renormalization-group theory based on model F starting with the calculation of the 
Green's function in Hartree approximation. We calculate the specific heat Cn(T, Q.), the correlation 
length ^(T, il), and the thermal-resistivity tensor pt{T,Q.) as functions of the temperature T for 
fixed values of the rotation frequency Q.. For nonzero f2 we find that all physical quantities are 
smooth near T\ so that the superfluid transition is a smooth crossover. We define a frequency- 
dependent transition temperature Ty, (fi) by the maximum of the specific heat and predict the shift 
Tx{^)-Tx = -1.2 Tx{2mAil^/TiY^'^" ■ For T < Tx(0.) we find mutual friction between the superfiuid 
and the normal-fluid component caused implicitly by the motion of vortex lines and calculate the 
Vinen coefficients B and B' . 
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I. INTRODUCTION 

Liquid *He in a uniformly rotating container is influ- 
enced by the rotation frequency $7 in two ways. First of 
all, a rotational flow is created which implies quantized 
vortices in the superfluid state for temperatures below 
T\ . Secondly, the centrifugal forces of the rotation cause 
a spatially dependent pressure P(r) which implies an in- 
homogeneity in the system. While usually the second 
influence is negligible, the rotational flow and the occur- 
rence of vortices plays an essential roleEl. 

In this paper we consider the critical behavior of uni- 
formly rotating liquid "^He in d = 3 dimensions for tem- 
peratures near Tx and present a renormalization-group 
theory, based on model F of Halperin, Hohenberg and 
Siggiaa. We will show that the rotation frequency fl is 
an external perturbation which drives the system away 
from criticality. According to FeynmanEI in uniformly 
rotating superfluid ''He straight vortex lines are present 
which are parallel to the rotation axis and which are uni- 
formly distributed in the helium. The circulation around 
a vortex line is quantized by 2'KTi/mi where is the 
mass of a ^He atom. Consequently, the density of vor- 
tex lines L with respect to the area perpendicular to the 
rotation axis is directly related to the rotation frequency 
by L = {m4/TTh)n. The absolute value L — {m4/TTh)Q is 
the total length of the vortex lines in a unit volume and 
hence is a measure of how many vortices are present in 
the helium. 

Close to the superfluid transition the critical fluctua- 
tions are very large. The correlation length ^ increases 
strongly if the temperature T approaches T\. However, 
for nonzero Q the mean distance between the vortex lines 
L^^/^ — [-Kh/m4^Y/'^ is another characteristic length of 
the system. We will show that the correlation length 
^ is bounded by this characteristic length according to 
^ < I/~'/^. Consequently, the thermodynamic quanti- 



ties like the specific heat are smooth functions of tem- 
perature near Tx. As a remnant of the critical singu- 
larities the quantities exhibit a maximum or an inflec- 
tion point at a temperature T\{yi) which is located below 
Tx. We will find a phase diagram for uniformly rotating 
''He which is shown qualitatively in Fig. |l|. The dashed 
line represents the transition temperature Tx{^) between 
the superfluid and the normal-fluid state. This temper- 
ature is not sharply defined, because for nonzero Q, the 
thermodynamic quantities are nonsingular. The tran- 
sition is smooth and shifted to lower temperatures by 
/^Tx{Vl) = Tx{fl)-Tx. Nevertheless, {T,fl) = (Ta,0) is 
the critical point in the phase diagram, for which the cor- 
relation length ^ diverges and the thermodynamic quan- 
tities are singular. 
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FIG. 1. The phase diagram of uniformly rotating liquid 
*He. The dashed line represents the transition temperature 
Tx{^) which separates the superfluid and the normal-fluid 
phase. The critical point is located at (T, Q) — {Tx, 0). 
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The macroscopic quantum coherence of the superfiuid 
component is described by the complex order-parameter 
field i)- For T > T\{Q) in the normal-fluid phase 
the order-parameter field '0 fluctuates around the average 
value (-0) = 0. On the other hand, for T < Tx{n) in 
the superfluid phase the order parameter = rj e^^ is 
nonzero. The amplitude rj fluctuates around a nonzero 
average value. However, the vortex lines in the uniformly 
rotating system imply a strong spatial variation of the 
phase if over all values between and 27r. Since the 
vortices move due to fluctuations, the phase Lp fluctuates 
strongly so that eventually the average order parameter 
is (■0) = also in the superfluid state. 

Conventional fleld-theoretic methods are not applica- 
ble in the present case because in the superfluid state 
symmetry breaking with a nonzero average order pa- 
rameter (0) is required for the construction of the per- 
turbation-, theory. However, recently we developed an 
approacha which can handle {tp) = also in the super- 
fluid state and which includes the effects of vortices in an 
indirect way. This approach starts with the calculation 
of the order-parameter Green's function in self-consistent 
Hartree approximation and includes the effects of criti- 
cal fluctuations by renormalization and application of the 
renormalization-group (RG) theory. In Sec. ^ we apply 
the approach to uniformly rotating *He. We describe the 
necessary modifications to include the rotation with fre- 
quency rJ. Since we consider also heat-transport phenom- 
ena, we additionally include an infinitesimal heat current 
Q as an external perturbation of the thermal equilibrium. 
Once the approach is set up, we can calculate several 
physical quantities explicitly. 

In Sec. m we calculate the entropy ^(r, Vl) and the 
specific heat Ca(T, f2) at constant rotation frequency f2. 
For nonzero we find that S and Co are smooth func- 
tions of the temperature T, so that the superfluid tran- 
sition is a smooth crossover. From the maximum of the 
speciflc heat we obtain the transition temperature T\{^) 
and calculate the shift ATA(f7) = Tx{^) ~ Tx- We con- 
flrm the phase diagram of uniformly rotating liquid '^He 
shown in Fig. |l|. The correlation length ^ = ^(T, fi) is 
calculated in Sec. We show that ^ is bounded by the 
mean distance between the vortex lines L~^/^ and flnd a 
maximum at T « Ta (fi) as expected. 

In Sec. ^ we consider heat transport in uniformly ro- 
tating liquid *He and calculate the thermal-resistivity 
tensor pt{T,^). For T < Tx{il.) in the superfluid state 
the thermal resistivity is strongly anisotropic and de- 
pends on the direction of the heat current related to the 
rotation axis. In Sec. VI we show that the thermal resis- 



tivity is caused by the mutual friction of the superfluid 
and the normal-fluid component due to the motion of 
the vortex lines, which was first observed by Hall and 
VinenH. We calculate the Vinen coefficients B and B' 
and compare the results with a previous theory for the 
motion of vortex linescl. Finally, in Sec. VII we compare 



the thermal resistivity of uniformly rotating superfluid 
^He with the thermal resistivity of nonrotating super- 



fluid "^He in the presence of a nonzero heat current Q, 
which was calculated previously in Rcf. ^ In this way we 
obtain the vortex density L of the turbulent superfluid 
flow induced by a flnite heat cujjxent Q which causes the 
Gorter-Mellink mutual frictioroa. 



II. RENORMALIZATION-GROUP 
THEORY FOR ROTATING *HE 

A. The model 

Dynamic critical phenomena in liquid ^He close to Tx 
are well described by model F which is given by the 
Langevin equations for the order parameter ■0(r, i) and 
the entropy variable m(r,t): 



dip 



SH 



SH 



dt " Sip* '^^^ 5m 



dt 
where 



dm . —nSH „ ^ /^SH\ 



(2.1) 
(2.2) 



H = J d'^r[iro(r)|0|2 + i|(V-zk)^|2-f^io|V'l' 

+ ^Xo^rn^ +lom\iP\^ - horn] (2.3) 

is the free energy functional and 9^ and 9m are Gaus- 
sian stochastic forces which incorporate the fluctuations. 
Because of the viscosity the normal-fluid compone nt is 
rotating uniformly with velocity Vn = $7 x r. In (2.3) 
the uniform rotation is incorporated by the wave vector 
k which is related to the normal fluid velocity by 



{m,i/Ti) Vn = [mi/h] {n X r) 



(2.4) 



We assume that the origin of the coordinate system is 
located on the rotation axis. The gravity and the cen- 
trifugal forces of the rotation imply a variation of the 
pressure in the helium according to 



P = Pa-p[gz-\{ny.TY] 



(2.5) 



where p is the density of liquid ^Hc and g = 981 cm/s^ is 
the gravitational acceleration. Since the critical temper- 
ature Tx — Tx{P) depends on the pressure P, it depends 
on the space coordinate r according to 

Tx{v) = Txo - (dTx/dP) p[gz - ^{n x r)^] (2.6) 

where! ^ = pg{-dTx/dP) = +1.273 /iK/cm. In the 
energy functional ( | 2.3|) the spatially dependent criti- 
cal temperature ( |2.6| ) is represented by the parameter 
To (r) = 2xo7o [^1 — Tx (r)] , where Ti is an arbitrary con- 
stant reference temperature. 

Thus, the rotation influences the ^He in two differ- 
ent ways. First, in the kinetic term of the energy func- 
tional (p.3[) the wave vector (2.4) related to the normal 
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fluid velocity generates straight vortex lines in the su- 
perfluid component. Secondly, the centrifugal force s im - 
ply a slightly space dependent critical temperature (2.6). 
In this paper we consider only the first kind of influ- 
ence which is first order in and which is the direct 
and most natural influence of the rotation. The second 
kind of influence, which is indirect via the pressure vari- 
ation and second order in Q, will be neglected. Further- 
more, we neglect gravity. In Sec. [II we will show that 
these neglections are justified for realistic experiments 
performed in a microgravity environment in space. The 
experiments with rotatims-superfluid ^He to determine 
the Vinen coefficientsEHl3"ll3 were performed with a slow 
rotation rotation up to two turns per second. For this 
reason we assume that the rotation frequency is about 
r2 « 27r s-\ 

For the calculation of the entropy and the speciflc heat 
in Sec. [II and of the correlation length in Sec. |^we may 
consider the system in thermal equilibrium, where the 
temperature T is constant. However, in Sec. ^ we con- 
sider heat transport and calculate the thermal-resistivity 
tensor in linear response. For this reason we must more 
generally consider the system in nonequilibrium, where 
an infinitesimal heat current Q is present which implies 
an infinitesimal temperature gradient VT. Thus, the 
temperature T(r) depends weakly on the space coordi- 
nate r. 

Model F is treated usually by field-theoretic means. 
The field-theoretic perturbation theory and the renornia.- 
lization-group (RG) theory were developed by Dohmtj. 
Our approach starts with the calculation of the order- 
parameter Green's function in Hartree approximation 
and is described in detail in our previous paper (the sec- 
ond paper of Ref . ^) . While the approach was originally 
developed for liquid "^He in the presence of a heat cur- 
rent Q, it is more generally valid and can be applied also 
to rotating helium with a few modifications. The main 
task is the explicit calculation of the Green's function, 
the set-up of the self-consistent equations for the effec- 
tive parameters, and the application of the RG theory. 



B. Unrenormalized Green's function in Hartree 
approximation and self-consistent equations 
for the effective parameters 

In Hartree approximation the self energy includes only 
the self-consistent tadpole diagrams and hence is inde- 
pendent of momentum and frequency. Consequently, 
the Green's function has the same structure as the free 
Green's function where only the parameters are replaced 
by the effective parameters. For our calculations we need 
the equal times Green's function 



G(r, 



4r' 



G(r,r') = (VXr,i)^*(r',i)) 



(2.7) 



In Ref. ^ we represented this Green's function in Hartree 
approximation in terms of the integral 



da e 



aA nB 



5[v- 



(2.; 



where 



A = -{Foh - (V - ik)2] - i(.go/2xo7o)Ar-o} , (2.9) 
B = -{nin - (V - ikf] + z(.9o/2xo7o)Aro} (2.10) 



are operators which are derived from the model- equa- 
tio ns (^T| )-( p.3[ ). Here k is the wave vector defined 
in (2.4) which incorporates the uniform rotation of the 
normal-fluid component, ri and Avq are the effective pa- 
rameters. In Ref. ^ it was shown that ri is related to the 
temperature parameter tq by 



where 



G(r,r) 



(2.11) 



(2.12) 



The parameters rg and Aro are related to the tempera- 
ture T(r) and the critical temperature Tx by 



ro = roc + 2x070 [^(r) - Tx]/Tx 
Aro = 2X070 [^(r) - To]/Tx 



(2.13) 
(2.14) 



where roc = in Hartree approximation and To is anarbi- 
trary constant reference temperature. (Compare (2.11)- 
( pl% with (3.31), (3.6), (3.30), and (3.26) in the sec- 
ond paper of Ref. ^, respectively). If r(r) and Tx are 
known, then ri, ro, andAro are determined by the self- 
consistent equations (2.11)-(2.14). Since we neglect the 
pressure variations by gravity and by centrifugal forces, 
Tx is constant. For the calculations of the entropy, the 
specific heat, and the correlation length in Sees. HI and 
IV the system is assumed to be in thermal equilibrium so 
that the temperature T(r) = T is constant, too. Thus, 
in this case the effective parameters ri, ro, and Aro are 
constant in space. 

For the calculation of the thermal-resistivity tensor in 
Sec. the presence of an infinitesimal heat current Q 
is needed. Consequently, there will be an infinitesimal 
temperature gradient VT which implies gradients of the 
effective parameters Vri, Vro, and V(A ro). S ince Tx is 
assumed to be constant, Eqs. ( ^.13[ ) and ( ^.14 ) imply 



Vro = V(Aro) = 2xo7o VT/Ta 



(2.15) 



On the other hand, from ( 2.11 ) we obtain 

Vri = Vro + 4uoVns . (2.16) 

Furthermore, we need a relation between the heat current 
and the gradients. For this purpose we take the average 
of the entropy equation (2.2) and obtain dt{m) + Vq = 
where 



\ dm I 



(2.17) 
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is the entropy current which is related to the heat current 
Q in physical units by q = Q/fceTA. The parameter Aro 
and the superfluid current Jg are defined byu 



Aro = 2x070 



(2.18) 



. om I 

J, = (Im[V'*(V-ik)?/;]) = 

= [(2^)-M(V-^k)-(V' + ^k')]G(r,r')],,^r ' (2-19) 

respectively. (Note that here the superfluid current is 
define d with respect to the rotating frame.) Thus, from 
(2.17) we obtain 



Q 



goksTx 



90 



Ao V(Aro) 

90 2xo7o 



(2.20) 



If the heat current Q is given, then the gradients V(Aro), 
Vro, Vri, and the temperature gradie nt V T are d eter- 
mined by the self-consistent equations ( ^.15 ), (2.16), and 



C. Renormalization 

The quantities rig, Vrig, and Jg, which appear in th e 
first order terms of the self-consistent equations (2.11), 
(2.16), and (2.20), exhibit infrared divergencies at crit- 
icality where ri ^ and — > 0. For this reason, the 
self-consistent equations must be renormalized and the 
RG theory must be applied to achieve a resummation of 
the infrared divergences and a proper treatment of the 
critical fiuctuations. We perform the renormalization in 
the same way as in Ref. ^ and use the concept of renor- 
malization by minimal subtraction of dimensional poles, 
which is described for model F in Ref. ^ The calcula- 
tions are performed at fixed dimension d — 4 — e (i.e. no 



e expansion)|_For the renormalization of (2.11) we need 
the relationsE^ 



'"o - roc = Zr-r , 



(2.21) 
(2.22) 



In Hartree approximation it is rpc = and the Z factors 
are given byo 



Zu = 1/[1 - 8u/e] 



(2.23) 



and Z^ = 1. Since ri is not renormalizedo, from ( ^.ll[ ) 
we obtain 



r = ri[l- 8u/e] - 4u{nyAd)ns 



(2.24) 



In the Appendix the Green's function G{r,r') is evalu- 
ated for infinitesimal gradients of the e ffectiv e parame- 
ters. Inserting this Green's function into ( ^.12 ) we obtain 



where J>(C) is defined by the integral 

poo 

M0= I 



(2.26) 



and £ — {%/ 2771/0)^/"^ is the characteristic length related 
to the mean distance between the vortices L~ ^J'^ . (The 
vortex density is i = l/27r£^.) Now, inserting ( 2.25 ) into 
(|2.24D we obtain 



ri£2{l + 8uA} 



(2.27) 



where 



A = 



n-l + e/2) 



{rieY'^^i+e/2irii^)-l ■ (2.28) 



Eq. ( 2.27 ) is the renormalized counterpart of ( 2.11 ). We 
have multiplied both sides by i'^ because ri'^ and ril'^ are 
dimensionless quantities. 



Next, by applying V to (2.27) we obtain the renormal- 
ized counterpart of ( 2.16[ ). We find 



Vri^ = {Vrif){l + 8uAi} 



where 



Ai = 



1 r 



(2.29) 



(2.30) 



The factors are multiplied on both sides of ( ^.29 ) be- 
cause Vr£^ and Vri^^ are dimensionless quantities. 



FOjLthe renormalization of (2.20) we need the rela- 
tionalj 

X070 = 7 {xoZmy^^Zri^iyAdy^^ , (2.31) 

5o-g(xo^™)'/'(A*7^<i)^/^ (2.32) 
Ao = XoZ^'X , (2.33) 

and furthermore Aro — ZrAr. The entropy current is 



renormalized by q = {xoZn 



11/2, 



Consequently, the 



Z factors cancel in the ratio q /gp — (qron/.9)(M'^/^d)^^^, 
so that the left-hand side of (2.20) needs not be renor- 
malized. On the right-hand side of ( 2.20| ) part of the Z 
factors cancel. Thus, we obtain 



Q 



Ad A V(Ar) 



Ai' 9 27 



{ZmZx 



(2.34) 



The dynamic couplings of model F like A and g oc- 
cur only in dimensionles s com binations like X/g. For 
this reason, we expresS|-(E^34r) entirely in terms of the 
dimensionless couplingal3 w = F/A, F = g/X, and 
/ = g^/XV = F^/wL In Hartree approximation the 
Z factors are given byQ 



Z„ZA = l/[l-//2e] . 
Thus, we rewrite ( ^.34| ) in the form 



(2.35) 



4 



Ad V(Ar) 



2e 



(2.36) 



The superfluid current Jg is obtained from ( 2.ig| ) by in- 
serting the Green's function G'(r,r'), which we have cal- 
culated explicitly in the Appendix. We find 



1 

7""r(-l + e/2 
F 



{M./2(r-i^') [e. X Vn] 



47W' 



^e/2('^i^')e.(e,. V(Ar)) 
-M^i^iT^t) [e, X [e, X V(Ar)]] 
+^;;2(^i^')[e.x V(Ar)]]} (2.37) 



where ^p(C) is defined in ( |2.26| ) and >lp(C), AAp(C), and 
A/'p'(C) are given by the integrals 



dv ■ 



^ - 1 j e-< 



shu Vthw 

2 



f + (u;"/w;')^ 7o 



dv 



sh^w 



cht) — cosf — -v\ 
\w' / . 

dv 



1 + («;"/«;')' h sh\ 
- w 



X 



— shu — sm — -v 

V \ w / 



(2.38) 



(2.39) 



(2.40) 



respectively. Now, inserting (2.37) into ( ^.36) an d re- 
placing Vri in terms of V(Ar) by using (|2"!29 ) and 
Vr = V(Ar), we eventually obtain 



e,(e, • V(Ar))/^3 
[e^ X [e^ X V(Ar)]]/^3 
V(Ar)]^3| (2.41) 
where Ai is defined in ( 2.30|) and 



1 + Ia' 

-{A"[e., X 



A' = - 
e 



"AC%(n£^) 



(2.42) 



e r(-l + e/2) 



47u>7F 
' 1 -h 8itAi 



A^,/2(ri£2) . (2.43) 



Final ly, w e renormalize ( p.l3 )-(2.15). Inserting ( p. 21 ) 
and (2.31) into (2.15) we obtain 



r = 2^{x^Z„:)'l\ii^lAafl\T(v) - T>]/Tx . (2.44) 

Clearly, in this equation the Z factor (xo-^m)"'^^^ does not 
cancel. For convenience we define the parameter 



27 



(2.45) 



so that (2.44) is written in the simple form 

r/^^ ^r-i[r(r) -TaI/Ta . (2.46) 



The renormalized counterparts of (2.14) and (2.15) are 
obtained analogously. We find 



Ar^2^r-i[T(r)-ro]/TA 



(2.47) 



^ v(Ar)//? = T-\^iTxy^VT , (2.48) 

respectively. 

D. Application of the RG theory 

By the renormalization a characteristic length scale is 
introduced which is described by the parameter fi. The 
RG theory is based on the fact that this length scale 
is arbitrary and may be changed according to fj, ^ ^il, 
where I is the RG flow parameter. As a consequence, 
the renormalized coupling parameters u{l), 7(0 j 
F{1), and /(/) depend on I. Furthermore, also the Z 
factors depend on /. (Note that the RG flow parameter 
I must be distinguished from the characteristic length 
i = (h/2 mifiy /^.) Now, the dimensionless parameter 
defined in (2.45) reads 



Adi^llY\^/^ 1 



\YnZ„Al)) 



XoZm{l)J 27(0 



(2.49) 



For convenience we will use r as the RG fiow parameter 
instead of I because r is closely related to the reduced 
temperature by (2.46) and the renormalized coupling pa- 
rameters u[t], 7[tJ, w[t], F\t], and /[t] were determined 
as functions of r in Ref. O. We identify iil — by 
the correlation length ^ = which in the asymp- 

totic region is given by ^(t) = ^0''" The identification 
/iZ — is corrjsct in one- loop order, corrections appear 
in higher orderaij. 

Now, we write the self-consistent equations for the ef- 
fective parameters in a form which is appropriate for the 
numerical evaluation. From (2.27), ( ^.29 ), and (2.41) we 
obtain 



ri£^{l + 8u[T]A} , 
(Vri£3){l + 8u[T]Ai} , 
Ad 



oe^ ^ 

gokBTx 2j[t]F[t] 



(2.50) 
(2.51) 

{[l + i/[r]^i]e,(e,.V(Ar))^3 



-[l + i/[r]A'][e, X [e,x V(Ar)]]^3 
if[T]A"[e.xV{Ar)]e} ■ (2.52) 



We replace fii nU = iji in ( [12^ ), ( ^!30|) , ( jli^ ), and 
( 2.43 ) and obtain the amplitudes 



5 



e 



WO' 



A' 



A" = - 



Lr(- 


-l + e/2) 




-WO' 


ir(- 


-l + e/2) 




-WO' 


ir(- 


-l + e/2) 




WO' 


r(- 


-l + e/2) 



KMnf ) - 1 



4:-/[t]w'[t]/F[t] 
1 + 8u[t]Ai 



,(2.53) 
(2.54) 
(2.55) 

(2.56) 



The functions J^p(C), MilC), A/"/,(0, an d AA" (0 are de- 
fined by the integrals (b.26) an d (|2.3§| )-( 2.40| ) as before, 
where however in ( ^.39 ) and ( |2.40| ) the parameters w' 
and w " m us t be replaced by w' [t] and w" [t] . Finally, 
from ( ^.4(: )-( 2.4^ ) we obtain the equations 



re 

Are 
Vre 



r-MT(r) 
r-i[T(r) 
V(Ar) e 



Tx]/Tx , (2.57) 
To]/Tx, (2.58) 
= T-\C/Tx)VT (2.59) 



which relate the effective parameters to the temperature 
T(r) and its gradient. 

Supposed the temperature T, the heat current Q, and 
the rotation frequency are known quantities, then Eqs. 
(2. 50)- (2. 59) are ten equations for the eleven unknown 
variables ri, r, Ar, Vri, Vr, V(Ar), A, Ai, A', A", 
and T. Thus, one variable remains undetermined which 
actually is the RG flow parameter r. In the spirit of the 
RG theory we must choose the flow parameter r so that 
in the perturbation series the infrared divergencies are 
resummed in an optimum way. By experience we find 
that the condition 



+ {32u[r]/e){l^e/2){^/er 



[nf + l]-<^/2 ^ I (2.60) 



is a good choice for fix ing r. 

The integrals ( |2.26D and ( |2.38| )-( |2.40| ) are well defined 
for C = rii"^ in the interval — 1 < C < +oo. For ( 
and for C, — 1 th e integrals can be evaluated asymptot- 
ically From (|^) we obtain 



T{p) 



1) 



-(p+1) 



for C > 
for C 



(2.61) 



Consequently, from ( 2.50| ) together with ( 2.53| ) we obtain 
riOil 



re 



{8u[T]/e)[{r,0)-'/'~m 

for ri^2 > +1 , 

-(16w[r]/e)(l - e/2){^/£f-'[riP + 1]-'/^ 

for ri^2 _^ „i . 

(2.62) 

For temperatures T well above T\ the correlation length 
^ decreases so that ^ <C ^. Eq. ( 2.60| ) imphes rii^ -t-1 



so that on the left-hand side of (2.60) the first term dom- 
inate s wh ile the second term can be neglect ed. T hus, 
Eq. ( 2.60 ) reduces into riO — 1 and Eq. ( 2.62 ) im- 
plies rr — rit"^. Consequently, for T ^ Tx the flow- 
parameter condition ([2.60|) reduces into 



re 



1 



(2.63) 



On the other hand, for temperatures T well below Tx it 
is ri£'^ —1, so that on the left-hand side of ( 2.60| ) now 
the second term dominates while the first term is small 
and negligible. Again it is ^ ^ i. Consequently, ( ^.62 ) 
implies that for T ^ Tx the flow parameter condition 
(2.60) reduces into 



2re = 1 



(2.64) 



Actually, ( pJ3| ) and ( pi are the standard flow-para- 
meter conditions of ReF|l| for T > Tx and T < Tx, 
respectively, where Q = and = 0. (Note that our 
r^^ is identified by r{l)/{nlY in Ref. O.) T hus, we have 



shown that our flow parameter condition ( 2.6C ), which 
is valid for finite and infinitesimal small Q, reduces 
to the standard flow parameter conditions of Ref. |l^ for 
temperatures T well above and well below Tx- For T near 
Tx it represents an interpolation. 

Now, Eqs. ( 2.50| )-(2.6C) are eleven equations which 
determine the eleven unknown variables uniquely. We 
solve these equations numerically for given temperature 
T, heat current Q, and rotation frequ ency f2 to determine 
the effective parameters. From (2.59) we obtain the tem- 
perature gradient VT. Once the effective parameters 
are known, physical quantities can be calculated explic- 
itly. This will be done in the next section. Since we 
consider liquid "^He in three dimensions, we set d = 3 
and e = 4 — d = 1 in all formulas when performing the 
numerical calculations. 

As an input we need the dimensionless renormalized 
couplings u[t\, 7[t], w[t] = w'[t\ -\-iw"[T], F[t], and /ix] 
as functions of r which have been determined by Dohmlij. 
Furthermore, we need the parameter which is related 
to the entropy at Tx- For liquid, helium at saturated va- 
por pressure this parameter isllj go = 2.164 x 10^^ s~^. 
To calculate the correlation length ^(r) = Cqt^" as a 
function of r we use the exponent v — 0.671 and the am- 
plitude ^0 — 1-45 X 10~* cm which were determined ex- 
perimentally in Refs. 15 and [l^. There are no adjustable 
parameters. 



III. ENTROPY AND SPECIFIC HEAT 

In our previous paperO the entropy was considered 
in Hartree approximation combined with the renorma- 
lization-group theory. The following result was obtained: 

S = Sx+t{B-V A[ (4i//a) + E[u*] ]r-" } (3.1) 

where 



6 



t = Tre = [T-Tx]/Tx 
E[u] = (2u)-i[l-Pi/p] • 



(3.2) 
(3.3) 



This formula is valid also for rotating "'He without al- 
teration. We just insert the effective parameters r and 
ri and the RG flow parameter r determined in Sec. ||. 
Since the entropy is an equilibrium quantity, the infinites- 
imal gradients may b e dis ca rded so that only the self- 
consistent equations ( |2.5C ), (2.53), and ( 2.57| ) must be 
solved. The specific heat Cn is then obtained by nu- 
merical differentiation with respect to the temperature 
according to 



'dS\ 



\dT/ n \ dt J n 



(3.4) 



where the rotation frequency Jl = ile^ is kept constant. 

The formula (3.1) is derived in the asymptotic regime 
close to Tx where ^ = ^qt'" and u[t] « u* = 0.0362. 
The critical exponents v = 0.671 and a — —0.013 and 
the nonuniversal amplitudes A = 2.22 J/molK and B = 
456 J/molK are obtained by compa ring t he t heoretical 
specific heat for f7 = defined, by (|3)-(i3 with the 
most recent experimental datall3, which were obtained in 
a microgravity environment in space. 




FIG. 2. The specific heat for rotating *He as a function of 
temperature. The solid line represents our theoretical result 
for Q = 2-K s"^ For comparison, the specific heat for £7 = 
is shown as dotted line. 

Realistic experiments with rotating ^He were per- 
formed for rotations about one turn per second (see Refs. 
1,01121). For this reason, we calculate the specific heat 
for Q^2Tr s'^. In Fig. | our numerical result is shown 
as solid line. For comparison the dotted line represents 
the specific heat at = 0. Clearly, for nonzero the 
specific heat is a smooth function of temperature with a 
maximum slightly below Tx. The position of the maxi- 
mum can be interpreted as the superfluid transition tem- 
perature Tx(0) for rotating ^He which, however, is not 



sharply defined. Thus, for nonzero Q the superfluid tran- 
sition is a smooth crossover. The shift of the transi- 
tion temperature ATx{^) — Tx{i^) — Tx is negative. For 
fl^2n s-i we find ATxi^) = -25 nK from Fig. |. For 
other rotation frequencies a scaling formula for ATx{^) 
can be derived. Scaling requires the relations = a 
and ATx{i^)/Tx = br where a and b are dimensionless 
constants. Inserting ^ — ^ot~'^ and £ = {h/2m4Qy/^ , 
eliminating t, and solving for ATx{^) we obtain 



ATxin) = ^MxTxi2mi^^n/hy/^'' 



(3.5) 



where Mx is a dimensionless constant of order unity re- 
lated to a and b. From the position of the maximum in 
Fig. ^ we find AIx — 1.2. The shift of the critical tem- 
perature (3.5) is represented in Fig. |^ by the dashed line. 
Thus, the phase diagram for uniformly rotating liquid 
""^He shown in Fig. ^ is confirmed. 

Our theory neglects the infiuence of the centrifugal 
forces which imply a sp atial variation of the critical tem- 
perature given by (2.6). For a sample of diameter D the 
maximum variation of the critical temperature is 



STx .max 

(n) = \dTx/dP\pn^D^ 



(3.6) 



To justify the neglection of the influence of the cen- 
trifugal forces we must require that (5T>,max(f^) is ten 
times smaller than ATx{fl) defined in (^?5|). Thus, for 
= 27r we req uire STx^maxi^) — 2.5 nK. Insert- 
ing this value into ( |3.6[ ) and solving for D we obtain 
the maximum diameter Umax of the sample for which 
the influences of the centrifugal forces may be neglected. 
For fl = 27r s ^^ w e obtain Umax — 0.6 cm, while 
from ( pTsj ) and (3_^) we find the frequency dependence 

realistic sample sizes. If the samples are made smaller so 
that D < Dmax is satisfied, then indeed the influence of 
the centrifugal forces may be neglected. 

On the other hand the influences of gravity on earth 
imply a much la rger variation STxg of the critical tem- 
perature. From (|2.6|) we obtain 



STxg = \dTx/dP\pgAz 



(3.7) 



where Az is the height of the sample. For Az « 0.6 cm we 
obtain STxg ~ 0.8 fiK which is much larger than the tem- 
perature scale in Fig. ^ For this reason, the maximum of 
the specific heat Cji induced by the vortices cannot be re- 
solved in an experiment on earth. The experiment must 
be performed in a microgravity environment in space. 

For temperatures well below Tx a first-order transition 
between a vortex liquid and a vortex lattice is expected, 
because for decreasing temperatures the fluctuations de- 
crease. However, our theory based on the Hartree ap- 
proximation cannot describe this transition. While the 
effects of vortices are included indirectly, our theory as- 
sumes the physical quantities to be homogeneous in space 
so that the system is always a vortex liquid. 
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IV. CORRELATION LENGTH 



where 



In the asymptotic regime close to T\{il) the correlation 
length is given by ^ = (a'O^^ ~ ^o''' Inserting the RG 
flow parameter t determined in Sec. ||, we calculate ^ = 
^(T, ri) as a function of temperature T for the rotation 
frequency fl — 2tt s^^. Our numerical result is shown 
in Fig. H as solid line. For comparison, the dotted line 
represents the correlation length at 17 = 0. Clearly, for 
nonzero Q the correlation length is a smooth function 
of temperature with a maximum slightly below Tx. This 
maximum is located close to the maximum of the specific 
heat in Fig. ||. 
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FIG. 3. The correlation length for rotating *He as a func- 
tion of temperature. The solid line is the correlation length 
{ for n = 27r s~^. The dashed line represents the magnetic 
length £. For comparison, the correlation length ^ for 0,-0 
is shown as dotted line. 

On the other hand, the characteristic length £ = 
{h/2m4fl)^^^ is shown in Fig. ^ as dashed line. By def- 
inition this length does not depend on the temperature. 
We find ^ < £ so that the correlation length is bounded 
by i. Near the maximum it is f ~ £. The character- 
istic length £ is closely related to the mean distance be- 
tween the vortex lines L~^/^. Because of L = l/27r€^ it is 
= (27r)-^/^^ « 2.5 £. Consequently, the correlation 
length ^ is bounded by the mean distance between the 
vortex lines according to ^ < L~^/^. 



V. THERMAL CONDUCTIVITY 
AND RESISTIVITY 



We eliminate the gradient V(Ar) from (2.52) and 
( 2.59 ) and solve the resulting equation for Q. Then we 
obtain 



Q - -AT4e.(e, • VT) 

+A;r(e, X (e, x VT)) 



goksAd {l + (/[r]/2)Ai} 



go fee 
goksAd 

td-2 



2j[t]F[t] 
{l + if[T]/2)A'} 

2j[r]F[r] 
inr]/2)A" 



2"f[T]F[T] 



(5.2) 
(5.3) 
(5.4) 



are three components of the thermal-conductivity ten- 
sor. While the dependence on the dimensio nalit y d is 
needed for the renormalization, the formulas (5.2)-(5.4) 
are evaluated for c? = 3. Clearly, Eq. (5.1) is a linear 



relation between the temperature gradient VT and the 
heat current Q. Solving this equation for VT we obtain 



Vr = -pT.ie,(e, -Q) 

+/9Jp(e^ X (e^ X Q)) - p'^{e^ x Q) 



(5.5) 



where pt,i, p't, and pi^ are three components of the 
thermal-resi stiv it y te nsor which are related to the con- 
ductivities (|l)-(|5l) by 



Pt 



Pt,i 
Vip'{ 



1/At4 
1/[A^ 4 



(5.6) 
(5.7) 



Clearly, Eqs. (qJ) and (^^) indicate that the heat cur- 
rent Q and the temperature gradient VT are not par- 
allel to each other. This fact is the consequence of the 
symmetry breaking by the rotation frequency = Oe^. 
In terms of the components the heat transport equation 



(5.1) can be written as 



gx + ^Q^/ = -(A^^+^A^^)(9,T^ 



idyT) 



Analogously, Eq. ( ^ ) can be written as 

dxT + idyT = -(pt + iPT)iQx + iQy 
dzT = -pT.iQz ■ 



(5.8) 
(5.9) 



(5.10) 
(5.11) 



\'i{B,^ X VT) (5.1) 



Thus, the heat transport decouples into two independent 
contributions, parallel and perpendicular to the rotation 
axis. The parallel heat transport (in the z direction) is 
very similar to the heat transport in nonrotating ^He: 
the thermal conductivity At,i defined in (5.2) has the 
same structure with the sa me a mplitu de A i as in the 
nonrotating case (compare (5.2) and ( 2^30 ) with (5.2) 
and (4.33) in the second paper of Ref.^, respectively). 
The only difference is the function Tp{C) which here is 
defined in a different way than in Ref. On the other 
hand, the perpendicular heat transport (in the xy plane) 
shows a completely different nature. In this cases the 
thermal conductivity AJ^ -I- iAJf and the thermal resistiv- 
ity yO^p -I- are complex. While the real parts AJp and 
pJp describe the dissipation of the heat current, the imag- 
inary parts \'J^ and p'J^ imply a coupling between the x 
and y direction so that the heat current and the related 
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temperature gradient do not have the same direction in 
the xy plane. The real part p'r^ may be interpreted as the 
longitudinal thermal resistivity, because for a given heat 
current in the xy plane it implies an antiparallel tempera- 
ture gradient. On the other hand, the imaginary part p'J^ 
may be interpreted as the transversal thermal resistivity, 
because it implies a temperature gradient perpendicular 
to the h eat cu rrent . The amplitudes A' and A" defined 
in ( ^.55| ) and ( 2.56 ) depend on compl etely differe nt inte- 
grals than Ai, which are defined in ( 2.38| )-(2.40) and in 
(2.26), respectively, and depend strongly on the rotation 
frequency f2. 

The special structures of the conductivities and resis- 
tivities reflect the fact that in rotating '^He there are 
straight vortex lines present which move and imply dis- 
sipation for a perpendicular heat flow but do not move 
for a parallel heat flow. The existence of a temperature 
gradient perpendicular to the heat current is the analogy 
of the Hall effect in electronic systems, where here the 
rotation frequency f2 plays the role of the magnetic field. 
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FIG. 4. The several thermal resistivities of rotating *He 
as functions of temperature for f2 = 27r s"'^. The longitudi- 
nal resistivity a^nd the transversal resistivity px for heat 
transport perpendicular to the rotation axis are shown as solid 
lines. The dashed line represents the resistivity Pt,i for heat 
transport parallel to the rotation axis. For comparison, the 
resistivity for f2 = is shown as dotted line. 



We calc ulate the thermal conductivities numerically by 
3 where the effective parameters and 



(|5^-(|5J) for d 

their gradients determined in Sec. ^ are inserted. From 

Pand ( pT?! ) we obtain the related resistivities. In Fig. 
r results for the thermal resistivities are shown for 
f2 = 27r s~^. The solid lines represent the longitudi- 
nal resistivity p'^ and the transversal resistivity p'-^ for 
the perpendicular heat flow, while the dashed line repre- 
sents the resistivity pT,i for the parallel heat flow. For 
comparison, the thermal resistivity for = is shown as 
dotted line. Clearly, for nonzero f7 the curves are smooth 
near T\. While /9t,i and p'^ increase continuously with 



increasing temperature, the transversal conductivity p',^ 
exhibits a maximum near T\. The inflection points of 
logj^qPx.i and logj^Q p'r^ and the maximum of logj^Q p'^ in 
Fig. ^ are located at a temperature below T\ close to the 
temperature of the maximum of the specific heat in Fig. 
1^. Thus, the shift ATx(il) of the superfiuid transition 
temperature for nonzero f7, defined in (3^), is observed 
also in the thermal conductivities where here the con- 
stant M\ is slightly different. 

For temperatures T well above and well below T\ 
asymptotic formulas can be derived for the thermal con- 
ductivities and resistivities. To do this we must evalu- 
ate the integrals J^p(C), 7Wp(C), AAp(C), and 7V^'(C) for 
C » -1-1 and for C — > —1, respectively. Wh ile for Tp{C,) 
the asymptotic formula is given by (2.61), analogous 
asymptotic formulas are found for the other integrals. 

First, we consider T T\ where C =LZi^^ ^ -1-1. In 
this case the fiow-parameter condition ( 2.60|) red uces to 
= 1. Then, for the amphtudes ( 2.53| )-( ^."5^ ) we ob- 
tain A « 0, Ai A' « -1/2, and A" w 0. Consequently, 
from (5.6) and (5.7) together with (^.2|)-( [5^ ) we obtain 
the thermal resistivities 



gokBAd {1 - /H/4} 



(5.12) 



and p'i^O where ^ ^ S.qt-" and t = (T - Tx)/Tx. Be- 
cause of (2.57) and (^.63|) the RG flow par amet er t is 



identifled by the reduced temperature. Eq. (5.12) is the 
well known formula for the thermal resistivity for f2 = 
and Q — > with the amplitude functions evaluated up 
to one-loop order (see Ref. |l^) . In Fig. ^ this resistivity 
is shown by the dotted line. Thus, for increasing tem- 
peratures T well above Tx the thermal resistivities px.i 
and p^p asymptotically approach the resistivity for f2 = 
where the imaginary part pi^ decreases and is negligibly 
small. This fact is clearly seen in Fig. ^. 

Secondly, we consider T <^ Tx where C — fi P — > — 1. 
In this case in the flow-parameter condition ( 2.60| ) the 
first term can be neglected so that 



r,e + 1 = {^(1 _ e/2){i/lf-^Y'^ (5.13) 



2/e 



which is very small for ^ <C ^. Because of ( 2.57 ) and 
(2.64) the RG flow parameter is related to the reduced 
temperature by r = 2{Tx — T)/T\. Inser ting t h e asy mp- 
totic formulas of th e i nteg rals ( |2.26| ) and (|2.38| )- (|2.40| ) for 
C ^ -1 into (|2^)-(|l5q) we obta in the amplitudes A, 
Ai, A\ and A". Then, from (0)-(^!4|) we obtain the 
thermal conductivities, and from (5.6D and (5.7) we ob- 
tain the related resistivities. Eventually, in leading order 
we obtain the asymptotic formulas 



PT,i ~ We){r,f + l)p'^ 
4^[t-]w'[t] 



Pt 



Pt 



gokBAd F[t] 



gokBAd F[t] 



Mr] {^/er 

8u[r] {i/ef 



(5.14) 
(5.15) 

(5.16) 
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where ^ ^ S.qt-" and I ^ {h/2miny/'^. 

For the heat transport perpendicular to the rotation 
axis for T <^ T\ the transversal and the longitudinal 
resistivities are related to each other by 



Pt/p't 



w"[t]/w'[t] 



(5.17) 



Neglecting the nonasymptotic effects due to the renor- 
malized coupling parameters 7[t], io'[r], i«"[t], and F[t], 
we obtain asymptotically in the leading order 



Pt Pt 



n (Ta - T) 



-l+a 



(5.18) 



where a = 2 — dh' = —0.013. Thus, in slowly rotat- 
ing superfluid ^He the thermal resistivities for the per- 
pendicular heat flow depend linearly on the rotation fre- 
quency. Since the resistivities are caused by vortex lines, 
this result is plausible: the resistivities must be propor- 
tional to the total length of the vortex lines per unit 
volume L which is related to the rotation frequency by 
L = l/27r^2 ^ (rn4/TTh)n. For n ^ the thermal resis- 
tivities p'rj^ and p'r^ vanish as expected. 

On the other hand, the heat flow parallel to the ro- 
tation axis is not directly influenced by the vortex lines, 
because the heat current flows parallel to the vortex lines. 
The dissipation of a parallel heat current and the related 
resistivity /9t,i must be caused by fluctuation effects. For 
this reason px,! is much smaller than the resistivity pJp 
of the perpendicular heat flow. This fact is clearly seen 
in F ig, [l and in Eq. ( 5.14 ) because the prefactor defined 
by (5.13) is very small. Neglecting the nonasymptotic 



effects due to the renormalized coupling parameters, we 
obtain asymptotically in leading order 



Pt,i 



-(d-2+i/e)v^-A/t 



(5.19) 



For d — i dimensions where e = 4 — d^lwe find /9t,i ~ 
^"^{Tx - TY'^", so that the resistivity for the parallel 
heat flow is quadratical in the rotation frequency. Again, 
Pt,i vanishes for 17 — > as expected. The temperature 
dependence of pT,i is very similar to the temperature 
dependence of the thermal resistivity in nonrotating "^He 
for nonzero heat currents Q (compare the dashed line 
in Fig. |4| with the solid lines in Figs. 3 and 4 in the 
second paper of Ref. ^). The related exponents for the 
approximate power laws for T <C are nearly the same. 



VI. MUTUAL FRICTION AND 
THE VINEN COEFFICIENTS 

Mutual friction in rotating superfluid "^He was inves- 
tigated first by Hall and VinenB. While originally the 
attenuatioa of second sound in resonant cavities was 
consideredH, the mutual friction implies dissipation and 
a related resistance for the heat transport in the super- 
fluid state. Here we show that the phenomenological the- 
ory of Hall and Vinen implies a heat-transport equation 



approach. In this way we show that mutual friction ac- 
cording to Hall and Vinen is derived from model F by 
our approach. 

Hall and Vinen proposed a mutual-friction force be- 
tween the superfluid and the normal-fluid component 
which is implied by the motion of the vortex lines in the 
presence of a superfiuid-normal-fluid counterflow. This 
mutual-friction force must be added to and subtracted 
from the two-fluid hydrodynamic equations for the super- 
fluid and the normal-fluid component, respectively. For 
a stationary homogeneous counterflow with the relative 
velocity Vg — v,i a relation for the temperature gradi- 
ent VT can be derived from the two-fluid hydrodynamic 
equations which reads 

(p/p„) sVT = -Bn {e, X (e^ x (v^ - v„)) 

+ (2-B')r!(e, X (v,-vn)) . (6.1) 

Here, B and B' are the Vinen coefficients which represent 
the phenomenological parameters of the mutual-friction 
force. The term with the coefficient 2 arises from the 
Coriolis force in the rotating frame. In the critical regime 
close to Tx it is p/pn ~ 1 and the entropy per mass is 
sKi sx = {h/m4){go/Tx)r^heTe go = 2.164 x 10" s'^ at 
saturated vapor pressurellJ. 

For a zero net mass current the superfluid-normal-fluid 
counterflow with relative velocity Vg — Vn is directly re- 
lated to the heat current Q. Considering the counterflow 
as a metastable state, in Ref. ^ the heat current Q was 
calculated for nonrotating superfluid "^He as a function of 
the counterflow wave vector k = {mi/h){vs — Vn). The 
calculation was performed in the critical regime near Tx 
using model F and the RG theory. Since the rotation fre- 
quency Q is very small, we may use this result here. For 
small counterflows and small heat currents the relation 
between Q and Vg — Vn is linear. While in Ref. |lj this 
relation was calculated in one-loop RG theory, for our 
purpose the relation in zero-loop RG theory is sufficient 
which reads 



Q 



Ad nu , 
8m H h ^""^ 



v„) 



(6.2) 



Now, we solve ( ^!^) for Vg — Vn and insert the resulting 
expression into (6.1). Then, we obtain the temperature 
gradient 



cd-2 

VT = I Su[t] Bn (e, X (e, x Q)) 
~'^KnAd 



g^KYiAd 
cd-2 



■8w[r] (2-B')r2(e2 X Q) . (6.3) 



g^k-B^ 

This formula has the same structure as ( |5.5| ). By com- 
parison we find the thermal resistivities 



which has the same form as (5.5) obtained within our 



(6.4) 
(6.5) 
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and pt,i = 0. For the heat flow perp end icular to the 
rotation axis the resistivities ( |6.4| ) and ( |6.5| ) are a direct 
consequence of the mutual friction due to the motion of 
the vortices. Because of = {nh/m4)L these resistivities 
are proportional to the vortex density L which precisely 
is the total length of the vortex lines per volume. Since 
p'rj^ implies a temperature gradient antiparallel to the heat 
current, the Vinen coefficient B describes the dissipative 
effects of the vortex lines. On the other hand, pJf im- 
plies a temperature gradient perpendicular to the heat 
current, so that the Vinen coefficient B' and the Coriolis- 
force coefficient 2 are related to reversible effects. For a 
heat flow parallel to the rotation axis the vortex lines do 
not move and hence do not cause dissipation, so that the 
resistivity px,! is zero. 

Within our approach which is based on the Hartree 
approximation combi ned with the RG theory we have 
derived the formula (5.5) for the temperature gradient. 
The terms involving the heat current perpendicular to 
the rotation axis have the same structure as the terms 



r'[r]/A[ 



'[r] = T"[t]/X[t], and F[t] = g[r]/X[T] 



on the right-hand side of (6.3). Thus, we conclude that 
our approach includes the mutual friction between the su- 
perfluid and the normal-fluid component caused by the 
motion of vortices. However, since the superfluid state 
is homogeneous in space and vortex lines are not treated 
explicitly, our approach includes the effect of the vortices 
indirectly. Furthermore, our approach includes also fluc- 
tuation effects, whi ch is reflected by the first term on the 
right-hand side of (^.5|) and the nonzero resistivity pT,i 
for the parallel heat fiow. This term is not present in 



Li 



We de termine the Vinen coefficients by solving (S.4) 
and (S.5) for B and 2 — B' and inserting the thermal 
resistivities p'j, and pi^ obtained within our approach. 
For small rotation frequencies Q the thermal resistivities 



are given by the asymptotic formulas (5.15) and (5.16). 
Thus, we obtain the Vinen coefficients 

B = {2mi/h)goreHr]w'[T]/F[T] , (6.6) 
2-B' = {2mi/h)g„Te47[T]w"[T]/F[T] , (6.7) 

which are defined in the limit 51 ^ 0. These formulas are 
well suited for numerical evaluations because all parame- 
ters and constants are known. The RG flow parameter is 
related to the reduced temperature by r = 2{T\ — T)/T\. 

Simpler expressions for the Vinen coefficients are ob- 
tained if we replace go by the renormalized counterpart 
g{l) — g[T\. From the renormalization equation (2.32) we 
obtain 



go^gil){xoZ^{l)Y'\{plY/AaY'^ 



(6.8) 



We solve (2.49) for (xo-Z^m(O)"'^^^ ^nd insert the resulting 
expression into ( |6.8| ). Then we obtain 

ffo=5M/(T^^27[T]) (6.9) 

where pi — has been identified. Now, we insert this 



result for g^ into (6.6) and (6.7). Furthermore, we in- 
sert the dimensionless renormalized couplings w'[t\ = 



Since most renormalized couplings cancel, we obtain 
eventually 



B = (4m4/?i)r'[T] 
B' = (4m4/?i)r"[r] 



(6.10) 
(6.11) 



This result is remarkably simple. Eqs. ( 6.1Cl| ) and ( |6.11 ) 
represent the final formulas for the Vinen coefficients ob- 
tained within our approach. 

The Vinen coefficients were calculated previously for 
model F within the renormalized mean-field theoryfl. In 
this previous approach the the vortex lines are consid- 
ered explicitly. The model-F equations are written first 
in the renormalized form and then solved as mean-field 
equations. Solutions were obtained which represent a 
single straight vortex line moving under the influence of 
a superfluid-normal-fluid counterflow. From the relation 
between the velocity of the vortex line vl and the relative 
velocity of the counterflow Vg — Vn the Vinen coefficients 
B and 2 — B' were extracted. Eventually, the effects of 
the critical fluctuations are included by application of 
the RG theory. The theoretical results were compared 
with the experimental data of Refs. |l0|-p^ obtained for 
temperatures Tx - T > 3 x 10"'' K. For B the agrefe 
ment between the experiments and the previous theoryQ 
is very good, while for 2 — B' there are some discrep- 
ancies. Consequently, the results of the previous theory 
may be viewed as correct and reliable. 

Close to criticality, the leading temperature depen- 
dence of the Vinen coefficients obtained from the previous 
theory is governed by 

B + i{2- B') « {2nn/h)T'[T]/[ci + C2h[T]F[T]f 

-ic3j[T]F[T]] (6.12) 

where ci = 0.2511, C2 = 0.391, and C3 = 0.795 (see Eq. 
(37) and Table 1 in Ref. ^, the winding number of the 
vortices is assumed to be n = 1). Since the product of 
the renormalized couplings 7[t]_F[t] is small, the term 
with the coefficient C2 can be neglected so that ( 3.12 ) 
simplifies into 



B 

2-B' 



(4TO4/7i)r'[T] X 1.99 , 
(4m4/;i)r'[T]7[r]F[r] x 6.30 



(6.13) 
(6.14) 



These formulas must be compared with ( 6.1(]| ) and ( |6.11 ). 
For B the leading critical temperature dependence is gov- 
erned by the renormalized coupling r'[r] in both cases, in 
the present theory and in the previous theoryQ. Quanti- 
tatively, the previous theory predicts a 1.99 times larger 
result for B than our present theory. On the other hand, 
for 2 — B' the leading critical temperature dependence is 
governed by different renormalized couplings for the two 
approaches, which is clearly seen in ( |6.11 ) and (6.14). 



The critical divergence of F" [r] is much weaker than the 
critical divergences of T'[t] and of r'[r]7[T]i^[T]. 

If we generalize model F by replacing the complex or- 
der parameter ip hy a complex vector ^E* = (^1, . . . ,tpn) 
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of n components, it turns out that the Hartree approx- 
imation is exact in the hmit n oo. The RG theory 
can be apphed but is not necessary in this case. Further- 
more, in the hmit — > for T < T\ the approximation 
reduces int o the mean-field theory. Consequently, Eqs. 
(6.10) and (3.11) may be interpreted as the Vinen coeffi- 
cients for n — oo \n renormalized mean-field theory. On 
the other hand, previouslyu we have calculated the Vinen 
coefficients for n = 1 in renormalized mean-field theory, 
which are given by (6.13) and (3.14) in leading order close 
to criticality. Thus, the Vinen coefficients allow a direct 
comparison of the physics for n = 1 and n = oo where 
the same approximation scheme is applied. While for B 
the agreement is quite good, for 2 — B' we find serious 
discrepances. 



VII. COMPARISON WITH THE GORTER- 
MELLINK MUTUAL FRICTION IN 
NONROTATING ^HE AT NONZERO 
HEAT CURRENTS 



lines. The therma l resi stivity ( [7.3[ ) has nearly the same 
structure as p'r^ in ( ^.15| ) where here the dephasing length 
^1 plays th e ro le of the characteristic length £. Thus, the 
resistivity ( |7.3| ) must be due to dissipation by the motion 
of vortex lines, which agree s wit h Vinen's picturdj. 

The thermal resistivities (5.14)-(5.16) represent the el- 
ements of the resistivity tensor for a heat flow in the 
presence of uniformly distributed straight vortex lines. 
For a tangle of vortex lines the thermal resistivity is ob- 
tained by taking the average over all directions of the 
vortex lines. Thus, we obtain 



p't 



PT,l) 



(7.4) 



which is just the trace of the resistivity tensor divided by 
d = 3. Because of pT.i <C resistivity pT,i of the 

heat fiow parall el to the vortex lines can be neglected. 
Now, inserting ( ^.15| ) into {7_A) we obtain 



Pt ~ t: 



2 tC 



d-2 



47 [t 



3 goksAd F[t] 



■ 8u[t] i^/e)' 



(7.5) 



Our approach was originally developed in Ref. ^ for 
nonrotating ^He in the presence of a nonzero heat current 
Q. For T < T\ we found mutual friction between the 
superfluid and the Bormal-fluid component according to 
Gorter and MellinkQ and calculated the Gorter-Mellink 
coefficient A. In the superfluid region we found a nonzero 
thermal resistivity which for temperatures T well below 
T\ is given by the asymptotic formulao 



Pt 



gokB F[t] 



( 8u[t] \^( Qe-' \ 



(7.1) 



where the RG flow parameter r is related to the reduced 
temperature by r = 2(Tx — T)/T\. Furthermore, we 
found a second correlation length ^1 which is a dephasing 
length for the order-paranietos held. For T ^ T\ we 
obtain the asymptotic formulaQ 



1 1/4 Ad gok-eTx 



8u[t] 



d-2 



(7.2) 



Here C is a dimensionless variable which depends loga- 
rithmically on the reduced temperature and varies be- 
tween —5 and —15 in the superfluid region. ( can be 
eliminated in favor of the dephasing length ^1. To do 
this, we solve ( |7.2| ) for ( and insert the resulting expres- 
sion into (7.1). Then, we obtain the thermal resistivity 



Pt 



re 



d-2 



47 [t 



gokBAd F[t] 



Mr] (2e/ei)' 



(7.3) 



Vinen arguedli that in superfluid "^He a nonzero heat 
current Q implies a turbulent superfluid flow and a tangle 
of vortex lines. The mutual friction and the dissipation of 
the heat current are caused by the motion of the vortex 



The prefactor 2/3 is due to the assumption that the direc- 
tion of the vortex lines is isotropically distributed. How- 
ever, since the isotropy is broken by a homogeneous heat 
current, the assumption may not be perfectly true. Thus, 
the correct prefactor will be slig htly different from 2/3. 

Now, we compare (7.5) with (7.3). We eliminate pt 
and solve the resulting equation for £^'^. Then, we obtain 
the total length of the vortex lines per volume 



L^l/i27rf)^{3/7r)C, 



(7.6) 



which is the density of the vortex lines induced by the 
heat current Q. This result indicates that the dephasing 
length ^1 can be interpreted as the mean distance be- 
tween the vortex lines. Inserting (7.2) for the dephasing 
length we find in leading order 



i - ~ (Vs - V„)2 



(7.7) 



so that the vortex density increases with the square of 
the counterflow velocity. Ka. (7.7) agrees with the vortex 
density proposed by VinenQ and with the ansatz of Gorter 
and MellinkQ for the mutual-friction force. Furthermore, 
in leading order we find that the temperature dependence 
of the vortex density is given by L ~ {Tx — T)2(<i-2)i'^ 



The average thermal resistivity (7.4) of the vortex tan- 
gle can be expressed in terms of the Vinen coefficie nt B 
and t he v ortex density L. For this purpose we insert ( |6.4[ ) 
into (7^), replace fl = (7r7i/m4)L, and obtain 



Pt = o 



2 ^ 



d-2 



3 g^ksAd 



r 1 7r?i „ 
8m T — BL 

m4 



(7.8) 



i.e. Pt ~ BL. Thus, there are two sources which deter- 
mine the magnitude of the thermal resistivity: the dissi- 
pation effects of a single vortex line described by B and 
the density of the vortices L. 
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In Ref . I we found that for a homogeneous heat current 
Q in superfiuid "^He our approach considerably overesti- 
mates the dissipation by the vortex lines. Our theory 
predicts a Gorter-Mellink coefficient A which iSjjiiiout 20 
times larger than observed in the experimentsDfl. This 
disagreement was confirmed recently by a measurementtS 
of the thermal resistivity px in superfiuid ^He close to Ta- 
We may now ask the question whether this large discrep- 
ancy arises due to the dissipation of the single vortex lines 
described by B or by the vortex density L. Rotating su- 
perfiuid ''He is the well suited system for a test of our 
approach to answer this question. Since the vortex den- 
sity L = {m4^/nh)Q is fixed and precisely known for this 
system, any observed discrepancy of the dissipation in ro- 
tating superfiuid ''He must be due to the Vinen coefficient 



B. In Sec. 



VI 



we have compared ourj)resent result for B 
with the result of the previous theoryB. While our present 
approach yields the correct critical temperature depen- 
dence, the values of B are smaller by a factor of 2 than 
expected. Thus, for a single vortex line our approach 
predicts a two times smaller dissipation than expected. 
On the other hand, for the heat transport at nonzero 
currents Q our approach predicts a 20 times larger dis- 
sipation. This large discrepancy can not be caused by 
the Vinen coefficient B. Rather, Eq. (7.8) implies that 
in Ref. ^ our approach predicts a 40 times larger vortex 
density L. Thus, we conclude that for superfiuid ^He in 
the presence of a heat current Q our theory considerably 
overestimates the density of the vortex lines L, while for 
the dissipative effect of a single vortex line described by 
B the correct order of magnitude is obtained. 

The discrepancy may possibly be explained in the fol- 
lowing way. The Hartree approximation implies that our 
approach effectively considers the generalized model F 
with n = oo complex order parameters. On the other 
hand, real liquid ''He is described by model F with n = 1 
complex order parameter. Thus, the discrepancy may be 
related to the difference of the physics between n = oo 
and n = 1. For n = 1 a homogeneous heat current Q 
and the related superfiuid-normal-fiuid counterfiow rep- 
resent a metastable state, which relaxes only by creation 
of vortices. For this purpose, energy barriers must be 
overcome, so that the vortex-creation rate, the vortex 
density L, and hence the dissipation are strongly sup- 
pressed. On the other hand for n = oo the heat current 
Q is unstable and vortices are created without an en- 
ergy barrier. Consequently, in this case the dissipation, 
the vortex-creation rate, and the vortex density L are not 
suppressed. Thus, the dissipation of the heat current and 
the vortex density may be considerably larger for n = oo 
(our theoretical approach) than for n — 1 (real superfiuid 
''He). The difference may be a large factor because the 
energy barriers occur in the argument of an exponential 
function. Thus, the quantitative discrepancy by a factor 
of about 20 or 40 can be explained in this way. 



VIII. CONCLUSIONS 

Our recently developed approach^, the Hartree approx- 
imation combined with the renormalization-group theory 
for model F, can be applied also to rotating ''He close to 
the superfiuid transition. For nonzero rotation frequen- 
cies n our theory predicts that all physical quantities are 
smooth and round near Tx- The superfiuid transition is a 
smooth crossover located at a temperature T\(r2) which 
is not sharply defined and shifted to lower temperatures 
by = Tx{n) - Tx = -MxTx{2m^eo^/K)^/^_- 

where Mx w 1.2. For tlitj-cptation frequencies « 27r s~^ 
of realistic experimentsOllHliS we find ATx(ri) « —25 nK 
which is very small. Thus, to observe the infiuence of the 
rotation on the superfiuid transition, experiments must 
be performed with a temperature resolution of a few nano 
Kelvins. Since on earth the gravity does not allow this 
temperature resolution these experiments must be per- 
formed in a microgravity environment in space. 

In superfiuid '^He for T < Tx the thermal resistiv- 
ity and the dissipation of the heat current strongly de- 
pend on the direction of the heat fiow related to the 
rotation axis. For a perpendicular heat fiow dissipa- 
tion is caused by the motion of the vortex lines in the 
superfiuid-normal-fiuid counterfiow. On the other hand, 
for a parallel heat fiow the thermal resistivity is consid- 
erably smaller and caused by fiuctuation effects. 

From the thermal resistivities of the perpendicular 
heat fiow we extract the Vinen coefficients B and B' in 
the limit — > 0. In this way we find that model F 
includes the mutual friction between the superfiuid and 
the normal-fiuid compoMnt, which originally was pro- 
posed by Hall and VinenU on phenomenological grounds. 
While our theory does not treat the vortices microscop- 
ically, we nevertheless obtain dissipation effects due to 
vortices. Consequently, our approach includes the vor- 
tices indirectly. 
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APPENDIX: CALCULATION OF 
THE GREEN'S FUNCTION 

The equa l time Green's fimction G(r,r') is defined by 



(2.8)-( 2.10 ). For the renormalization we need addition- 

-1 1/2 

ally the relations Tq = F and ip = Z^ 'i'ncw where, 
however, in Hartree approximation it isQ Zy — \ and 
Z^ = 1. Thus, from (^|^) we obtain the renormalized 
Green's function 



G(r,r') = 4r' 



da e"^ e"^ 8{v 



(Al) 
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By using the relations ( p. 31), (p.32 ), and Aro — ZrAr, 
in the operators ( p^ ) and ( 2.10| ) we replace the bare pa- 
rameters by the renormalized counterparts. All Z factors 
cancel except Z r wh ich, however, is unity here. Thus, 
from ( p^gj ) and ( 2.10 ) we obtain 



A = -{r[ri - (V - zk)2] - i{g/2^)Ar} , (A2) 
B = -{r*[ri - (V - ikf] + i{g/2'y)Ar} . (A3) 

In the following we must figure out, how the operators 
e""^ and e"^ act on the delta function 5{r—r') to evaluate 



the integrand in (Al) 



We derive an explicit integral formula for the Green's 
function for rotating "^He in the presence of an infinites- 
imal heat current Q and an infinitesimal temperature 
gradient VT. For the effective parameters ri and Ar we 
use the ansatz 



ri = ai + bir , 
Ar = a + b r , 



(A4) 
(A5) 



where ai, a are constants and bi, b are infinitesimal 
vectors which are related to the tempe ratu re gradient 
VT. We evaluate the Green's function (Al) up to first 
order in bi and b. For this purpose we decompose 



A = Ao 



AA 
AB 



where 
Bo 

and 



-{T[a^-[V -i\^f]-i[gl2^)a] , 
-{r*[ai-(V-zk)2]+i(g/27)a} 



A^ = -[rbi-i(g/27)b]r , 
AB = -[r*bi + i{g/2'y)h] r 



(A6) 
(A7) 



(A8) 
(A9) 



(AlO) 
(All) 



We expand the operator e"'^ 
order in AA according to 



1 



a 



n+l 



n— ^ ^ 



"■^''(aAA)e-^"-^''} 

[Ao,A.4]...]}e"^« 



\A. 



(A12) 



where the number of operators in the multiple com- 
mutator is n. Analogously, we expand e"^ = QO'Bo+aAB 
up to first order in AB according to 



aB 



s"^°|l + ^ dA e-^"^°(aAB)e^"^"} 



n+l 



\AB. B, 



oj 



■Bo]} 
(A13) 



where the number of operators Bq in the multiple com- 
mutator is n. 

Next, we evaluate the multiple commutators explicitly 
and resum the series' over n by the hyperbolic functions. 
For convenience we define the vectors 



S = bi - z(.g/27r)b, 
S* = bi+z(g/27r*)b 



(A14) 
(A15) 



so that AA = -FSr and AB = -F*S*r. The rotation 
frequency = fie^ is incorporated via the wave vector 
k, which is defined by (|]J) and which represents the uni- 
form rotation of the normal-fluid component. Eventually, 
we obtain the operators 

^ {1 _ aTSr - (aF)2(e^S)(e,(V - ik)) 



-i[ch(2Fa/£2) _ ^] ^4 ^ ^ s))(^ _ .j^) 
■|[sh(2Fa/£2)- (2Fa/<'2)]^4 (e^ x S)(V-ik)} 



(A16) 



and 



nB 



„aB, 



» {l - aF*S*r + (aF*)2(e^S*)(e^(V - ik)) 
-^[ch{2T*a/f) - 1] £^ {e, x {e, x S*))(V - ik) 
-i[sh(2F*a/£^) - {2T*a/e^)] £^ (e^ x S*)(V - ik)} 

(A17) 

where £ = (?i/2m4ri)^/^. Now, we insert thes e operators 
into the formula for the Green's function (Al) and obtain 



/>oo 

G(r, r') = 4F' / da {■ ■ ■} e"^" e"^° {• • •} S{r - r') . 
Jo 

(A18) 

where the first curved brackets {• • •} are iden tified by 
the expressions in the curved brackets in ( A16| ) and the 
second curved brackets {• • •} are ide ntifi ed by the ex- 
pressions in the curved brackets in ( A17 ). Because of 
[Ao,Bo] = we find 

gOAo gOBo ^ ^a{Ao + Bo) ^ g-a2r'[ai-(V-!k)^] ^ (A19) 



In ( A18 ) in the second curved brackets the operator V 
acts directly on the delta function so that we can use the 
identity 



V6{r 



-V'S{r 



(A20) 



Analogous identities are valid for r and k. Thus, we 
obtain the Green's function 



/•OO 

G(r,r') = 4F'/ da {■■■}{■■ ■}' 

Jo 

xe-"2r'[ai-(V-.k)^]^(j._^/ 



(A21) 



where in the second curved brackets {• • •}' the operators 
V, r, and k are replaced by — V, r', and k', respectively. 
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(k' is defined by (2.4) witli r replaced by r'.) We eval- 
uate the product {• • •}{• • •}' up to the terms first order 
in S and S*. Furthermore, we substitute 2T'a — > a in 



the integral. Then, from (A21) we obtain the Green's 
function 

G(r,r') = 2/ da (l - — Sr - — S*r' 

-(^)'(e.S)(e,(V-zk)) 
/ar*N2 



(^^j (e.S*)(e,(V' + zk')) 

V rt \ 1 

' £4(e, x(e, xS))(V-ik) 



ch 



r a 



e*ie, X (e, X S*))(V' + ik') 



M^)-Q)]^VxS*)(V' + zk')} 
X e-^I^i-C^-'k)'] 5(r - r') . (A22) 



Since k defined in (|^) has the same form as the vector 
potential of a homogeneous magnetic field, the eigenfunc- 
tions and eigenvalues of the operator (V — ik)^ are the 
Landau levels where i — (h/2m4^ny^^ is t he "magnetic 
length" . The factor in the last line of ( A22 ) can be eval- 
uated by representing the delta function in terms of the 
Landau-level eigenfunctions. Replacing the exponential 
operator by the related eigenvalues and performing the 
sum over all eigenvalues, we obtain eventually 



X expl —aax 



(47ra)'^/2 sh(a/^2) 

p2 aie e 



4ath(a/^2) 4q! 



i- (A23) 



where = {x — x'Y + (y — v'Y and = (r — r')^ — . 

In the next step we apply the dif feren tial operators V 
and V in the cur ved b racket s of (A22) onto the factor 
( |A23| ). We insert ( |A1^ ) and ( |aT5|) for S and S*. Since 
the dynamic parameters of model F occur only in di- 
mensionless ratios, from now on we write the formulas 
in terms of the dimensionless renormalized parametersli3 
w ~ T / X, w* = r*/A, and F — g/\. Then, after some 
manip ulations of the terms in the curved brackets in 
( A22 ) by using the addition theorems of the hyperbolic 
functions and by using F/F' = w/w' = 1 -|- iw" /w' and 
F*/F' = w* /w' = 1 ~ iw" /w', we obtain the following 
result for the Green's function: 

G(r,r') = e-'^=('-^'^')/2^' 



(4 



2 f°° da 



F 



+ 2(e,.(r-r'))a^(e,.b) 



P r/ 

+ *(r-r')y[( 



1 



F 



2w' 



{gz X bi) 



^sh(^)-sin(^)](e.xb))]} 



a/e 



sh{a/P) 



expl ~aai 



a/e 



4ath(a/^2) ■ 
(A24) 



Finally, we use the formula « {1 -f x} for small x to 
combine the terms in the curved brackets into an expo- 
nential function. Thus, eventually we obtain 



G(r,r') 

x exp| 
where 

and 



da ajP' 



(47r)'^/2 7o oi'^l'^ sh(a/^2) 



-ar\ 



4a th(a/^2) 



C_\ g»K(r-r') (^25) 

4a/ 



ai -)-bi(r-|-r')/2 



(A26) 



F 



K = a- e^(e2 • b) 

87W' 

aje 



P2 r 



th(a/^2) 
F 2w' 



(e^ X bi) 



sh(a/^2) 4^ w'2 +^t,//2 



cos 



(e^ X (e^ X b)) 



w , /a\ . (w a\i , , ,\i , . 



an d C are de fined below (A23). Eq. (A25) together 



with ( |A26| ) and (^2|) is the final result for the Green's 
function G(r, r'), which is valid for infinitesimal gradients 
bi — Vri and b = V(Ar) up to first order in bi and b. 
For K = and fi = ri = const. Eq. ( A25) represents the 
Green's function for rotating "^He in thermal equilibrium. 

The wave vector K defined in ( |A27 ) is the sum of the 
two terms of different natur e wh ich are orthogonal to 
each other. The first term in (A27) represents the contri- 
bution parallel to the rotation axis and does not depend 
on the rotation frequency 17. Th e sec ond term, which is 



given by the remaining terms in (A27), is in the xy plane 
perpendicular to the rotation axis and depends on via 
the characteristic length I. In the limit — > and £ ^ 00 
the wave vector reduces to K — > a(F/87w')b. This wave 
vector is well known from the theory of heat transport 
in nonrotating ^He (see e.g. (A28) in the second paper of 
Ref. a. 
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While the scalar and vector products involving are 
defined in d = 3 dimensions, our formula for G(r, r') is 
valid for arbitrary d dimensions. The scalar and vec- 
tor products in (A25) and (A27) must be replaced by 
appropriate generalized expressions. This, however, is 
straightforward: —e^ x (e^ x b) is replaced by the pro- 
jection of b onto the xy plane and ez{ez ■ b) is replaced 
by the projection of b onto the subspacc perpendicular 
to the xy plane. Furthermore, x bi and x b are 
replaced by the projections of bi and b onto the xy plane 
combined with a rotation of angle tt/2. 

The density rig and the superfluid current Jg ar e ob- 
taine d by inserting the Green's function ( A25| ) into ( ^.12 ) 
and ( 2.19| ), respectively. We find 



(4^)<i/2 



da a/i'^ 



sh(a/^2) 



and 



da ajP' 



(47r)'^/2 7o a'^/z sh(a/£2) 



where K is defined in (|A27D . We replace 



(47r)''/2 



-Ai— 

e r(-i + 



e/2) 



(A28) 



(A29) 



(A30) 



where is the geometric al fa ct or wh ich occurs in the 
renormalization equations ( |2.22| ), ( 2.31 ), and ( 2.32 ). Fur- 
thermore, we substit ute a j^^ v. T hus, eventually we 
obtain the formul as (^.25 ) and ( ^.37 ) with the integrals 
(p6|) and (|238|)-(|2^q)~"^ 
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